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ABSTRACT 

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift 
range z ~ 7 — 12 using deep images from the 2012 Hubble Ultra Deep Field (UDF12) campaign, data 
which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased 
signal-to-noise ratio offered by the UDF12 imaging to provide improved size measurements for known 
galaxies at z ~ 6.5 — 8 in the HUDF. Specifically, we stack the new deep F140W image with the 
existing F125W data in order to provide improved measurements of the half-light radii of z'-band 
dropouts (at z ~7). Similarly we stack this image with the new deep UDF12 F160W image to obtain 
new size measurements for a sample of F-band dropouts (at z c± 8). Secondly, because the UDF12 
data have allowed the construction of the first robust galaxy sample in the HUDF at z > 8, we 
have been able to extend the measurement of average galaxy size out to significantly higher redshift s. 
Restricting our size measurements to sources which are now detected at > 15a, we confirm earlier 
indications that the average half-light radii of z ~ 7 — 12 galaxies are extremely small, 0.3 — 0.4 kpc, 
comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm 
that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first 
evidence that this trend continues beyond z c± 8. Modeling the evolution of the average half-light 
radius as a power-law, oc (l + z) s , we obtain a best-fit index of s = —1.28 ±0.13 over the redshift range 
z ~ 4 — 12, mid- way between the physically expected evolution for baryons embedded in dark halos 
of constant mass (s = — 1) and constant velocity (s = — 1.5). A clear size-luminosity relation, such as 
that found at lower redshift, is also evident in both our z- and F-dropout sample. This relation can 
be interpreted in terms of a constant surface density of star formation over a range in luminosity of 
0.05 — 1.0L* =3 . Our results also strengthen previous claims that the star- formation surface density in 
dropout galaxies is broadly unchanged from z ~ 4 to z ~ 8 at Esfr — 2 M yr _1 kpc -2 . This value is 
2 — 3 orders of magnitude lower than that found in extreme starburst galaxies, but is very comparable 
to that seen today in the centers of normal disk galaxies. This provides further support for a steady 
smooth build-up of the stellar populations in galaxies in the young universe. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift — galaxies: 
structure 



1. INTRODUCTION 

Considerable progress has been made in charting the 
abundance of galaxies at z ~ 7 — 10 from deep imaging 
with ground-based observations and various campaigns 
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undertaken with the Wide Field Camera 3 infrared chan- 
nel (WFC3/IR) on Hubble Space Telescope (HST). Sam- 
ple selection makes use of the well-established dropout 
technique, which takes advantage of the unique spec- 
tral characteristics of high-redshift star-forming galax- 
ies, i.e., a blue UV spectrum and a sharp drop in flux 
at wavelengths shorter than Ly<x These complemen- 
tary studies have identified a large number of dropout 
galaxies at z ~ 7 and beyond. Investigating the abun- 
dance of dropout galaxies over 7 < z < 10 has re- 
vealed a clear decrease in the number density of luminous 
galaxies with increasing redshift (e.g. , lOuchi et al.ll2009l: 
McLure et all 120101: IC astella no et all 120101 : lOesch et al.l 
mObtlBouwens et al.ll20llbh . 

Characterizing the evolution of galaxy morphologies 
and sizes is useful for understanding galaxy formation 
history. Analytical studies have calculated the size- 
redshift relation of disk galaxies, suggesting the typi- 
cal size of galaxies of a given lu minosity is expected t o 
decrease with increasing redshift (|Mo et "all 1 19981 fl999h - 
The virial radius of a dark matter halo scales with red- 
shift and virial velocity or virial mass. Assuming that the 
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exponential scale length of the baryonic disc scales with 
the virial radius, the sizes of disks are expected to scale 
with redshift, proportional to H(z)~ 1 at a fixed mass or 
H(z) ~ 2 / 3 at a fixed circular velocity (e.g. JFerguson et al.l 
l2004f ). where H(z) is the Hubble parameter which scales 
as ~ (1 + zfl 2 at high redshifts. 

Earlier observations have reported that the sizes (half- 
light radii) of dropout galaxies decrease according to 
about (1 + z)- 1 up to z ~ 7 (jOesch et all l2Q10ah . 
which is expected at fixed halo masses, and in good 
agreement with prev ious estimates at lower redshifts 
(jBouwens et a l. 2004). However, because they used only 
the first epoch of their survey data with the WFC3/IR, 
their analysis still shows large uncertainties, especially in 
their fainter sample. Therefore, it is also consistent w ith 
(1 + z)- 1 ' 5 (jFerguson et al.l 120041 : lHathi et all 120081 ). as 
e xpected for s i zes tha t scale with halo circular velocity. 

lOesch et al.l (j201 0a) have also reported that the star- 
formation rate (SFR) surface densities of dropout galax- 
ies remains constant from z ~ 7 to z ~ 4. They suggest 
a possible explanation is that the average star formation 
efficiency is very similar in all these galaxies, and that 
feedback effects change the mode of star formation by 
only a small amount. It would be interesting to see if 
this possible trend continues toward higher redshifts, to 
infer star-formation activities in galaxies at earlier epochs 
of galaxy formation. 

Recently, a new campaign was carried out with the 
WFC3/IR to significantly deepen the Hubble Ultra 
Deep Fi eld in 2012 (GO 124 98; PI: R. Ellis, hereaf ter 
UDF12: lEUis et al.ll20l2l and IKoekemoer et al.ll2012l for 
the project description); this yields the deepest near- 
infrared images ever obtained. Additiona l scientific re- 
sults f rom this projec t are p resented in iDunlop et all 
(120121). ISchenker et al.l ([20 121 ). iMcLure et al.l (|2012l ). and 
iRobertson et al.l (|2012[ ). In this paper, we study mor- 
phologies of z ~ 7 — 12 galaxies based on the complete 
WFC3/IR UDF12 data set. The advantages of the new 
images are (i) a new F140W image and a deeper F160W 
data from which we obtain robust estimates on the rest- 
frame UV morphologies of galaxies not only at z ~ 7, but 
also z ~ 8 — 12 for the first time, (ii) a deeper F105W 
image which enables us to safely exclude contaminations 
by foreground sources from our galaxy samples at z ~ 8 
and beyoncQ. The purpose of this paper is to investigate 
the galaxy size and SFR surface density evolution beyond 
z ~ 7, and the correlation of size with UV luminosity. 

The outline of this paper is as follows. After describ- 
ing the imaging data used in this study in Section [2j we 
summarize our dropout galaxy samples in Section [3l Our 
size analysis is described in Section HJ In Section [5j we 
investigate the size-luminosity relation and the size evo- 
lution and discuss the implications. A summary is given 
in Section [6l Thro ughout this paper, we use magnitudes 
in the AB system (|Oke fc Gundfl983l ) and assume a flat 
universe with (O m , Qa, h) = (0.3, 0.7, 0.7). In this cos- 
mological model, an angular dimension of 1.0 arcsec cor- 
responds to a physical dimension of 5.365 kpc arcsec -1 
at z = 6.7, 4.818 kpc arcsec -1 at z = 8.0, 4.465 kpc 
arcsec -1 at z = 9.0, and 3.683 kpc arcsec -1 at z = 11.9. 

11 We do not use our deep F105W image for the morphology 
analysis of z ~ 7 galaxies, since a redshifted Lya and the continuum 
break of an object at z > 6.4 enters the F105W band. 



We express galaxy UV luminosities in units of the char- 
acteristic luminosity of z ~ 3 galaxies, LZ—*, whi ch cor- 
responds to Mieoo = -21.0 (jSteidel et al.lll999f ). The 
four WFC3/IR filters we use, F105W, F125W, F140W 
and F160W are denoted by Y105, J125? -^140 and i?i6o, re- 
spectively. We also use four ACS filters, F435W, F606W, 
F775W and F850LP, which are denoted by ,8435, V 606 , 
Z775 and ^850, respectively. 

2. OBSERVATIONS 

The primary data set used in this morphology analy- 
sis for z ~ 7 — 12 galaxies is the ultra-deep WFC3/IR 
observations taken for the UDF12 campaign combined 
with images taken for the UDF09 campaign (GO 11563; 
PI: G. Illingworth). In the UDF09 campaign, WFC3/IR 
data was obtained over three fields: the HUDF main, 
and two parallel fields. The UDF12 campaign has ob- 
tained 128 orbits of WFC3/IR data over the HUDF 
main field. We have combined all the exposures in- 
cluding the data from other HST programs (GO12060, 
12061, 12062; PI: S. Faber, H. Ferguson; GO12099; PI: 
A. Riess). In total, the observations over the HUDF main 
field include 253 orbits (F105W: 100 orbits, F125W: 39 
orbits, F140W: 30 orbits, F160W: 84 orbits). A more 
de tailed description o f the UDF12 data set is provided 
by IKoekemoer et al.l (|2Q12h , and the final reduced data 
are being made publicly available as High-Level Science 
Productq3 that are delivered to the Space Telescope Sci- 
ence Institute archive, and further details and current 
updates about the survey are provided at the project 
websit<£EJ 

To minimize the effects of morphological if-correction 
and take the advantage of the UDF12 campaign, we mea- 
sure sizes of galaxies in the images of the WFC3/IR band 
that is the closest to the rest-frame 1600 — 1700A. A stack 
of the PSF-matched J125- and Ji4o-band images is used 
for zgso-dropouts, a stack of the PSF-matched Juo- and 
i7i6o-band images for Yios-dropouts, and the i^i6o-band 
image for candidates at z > 8.5. Their 5<j limiting mag- 
nitudes are 29.8 (J125 + Juo), 29.7 (Juo + #160), and 
29.5 (Hiqq) within filter-matc hed apertures, w hich are 
0.45 - 0.50 arcsec in diameter ((Ellis et al.ll2012l ). We use 
images with a pixel scale of 0.03 arcsec pixel -1 . 

3. SAMPLES 

We investigate the sizes of z ~ 7 — 12 galaxies based 
on the the z ~ 7 — 8 samples selected by iSchenker et al.l 
(|2012[ ) and the z > 8.5 samples whose photomet- 
ric redshifts from SED fittin g analysis a re available in 
IMcLure et all (|2012f ) (see also lEllis et aLll2012[ ). Here we 
briefly summarize how these galaxies are selected. 

To select star-form ing galaxies at z ~ 7 — 8, 
ISchenker et al.l (|2012[ ) applied the dropout technique, 
which probes a blue UV spectrum and a spectral break 
blueward of Lya due to IGM absorption. For z ~ 7 zsbo- 
dropout galaxies, they first required a 3.5a detection in 
Y105 plus one of the other filters which probe longer wave- 
lengths (J125, Juo, or Hiqq). Then they applied the two 
color criteria: zs^o — Y105 > 0-7 and Y105 — J125 < 0-4. 
Also the following criteria were used; (i) the significance 

12 http: / /archive. stsci.edu /prepds /hudf 12/ 

13 http://udfl2.arizona.edu/ 
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Table 1 

Bright ^850" and Y^q^ -Dropout Galaxies in the HUDF Reported in the Literature 
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Note. — (1) Bouwens et al. (2008) ; (2) Oesch et al. (2009); (3) Oesch et al. (2010b) and Bouwens et al. (2010); (4) McLure et al. (2010); (5) 
Finkelstein et al. (2010); (6) Yan et al. (2010); (7) Bunker et al. (2010); (8) Wilkins et al. (2011) and Lorenzoni et al. (2011); (9) Bouwens et al. (2011b) 
and Bouwens et al. (2011a); (10) McLure et al. (2011). UDF 12-4258-6567 has been spectroscopically confirmed by [Font ana et al. l20ll]). 

t 1 Close to UDFz-39557176. 
t 2 Close to zD4. 



Table 2 

Photometric Redshifts for Bright zsso- an d Yios-Dropout Galaxies Found in the HUDF 
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-4470- 


-6442 
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7.32 
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■3952- 


-7173 


1422 
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(7.00-8.05) 


2055 
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(7.48-8.24) 
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■4314- 
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1678 
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7.30 


UDF12- 


■3722- 
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1574 


7.20 


(6.55-7.60) 
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7.40 


(7.12-7.70) 


40410 


7.70 


UDF12- 


■3813- 


-5540 


1721 


8.45 


(7.75-8.85) 


125 


8.61 


(8.19-8.89) 


70020 


8.30 



Note. (i) iMcLure et~aD (I2010D : (ii) Fin kelstein et al.l (I2010D : (iii) Mc Lure et"aD (120121 ) . 

5z p hoto denotes lcr confidence interval. 

^ 1 This object is not included in the catalog of McLure et al. ( 2012), because its photometric redshift 
is less than 6.5. 



is less than 2.0a in .8435, ^06, and ^775; (ii) the sig- 
nificance is not more than 1.5<j in more than one band 
among ,8435, Veoe, and ^775; (iii) Xo P t < 5.0 E3- For z - 8 
Yios-dropout galaxies, they required a 3.5a detection in 
J125 and one of the other filters which probe longer wave- 
lengths, J140 and #160- From the detected objects, they 
selected dropouts which satisfy the two color criteria, 
^105 — -^125 > 0.5 and J125 — #160 < 0.4, and the fol- 
lowing criteria for the optical data; (i) the significance is 
less than 2.0a in ^435, ^606, ^775, and (ii) the sig- 
nificance is not more than 1.5a in more than one band 
among £435, ^606, ^775, an d z 85 o; (hi) xl v t < 5.0. By us- 
ing these selection criteria, iSchenker et al.l (j2Q12h identi- 



L4 Xopt is defined by xl pt 



E,SGN(/i) Ui/GiY where is the 
?i is the uncertainty of fi, and SGN(/^) is 1 if 



flux in band 

fi > and —1 if fj < 0, considering the bands shorter than Lyo: 
(Bouwens et al. 2011b). For Z85o-dropouts, B435, V606, and ^775 
are considered, and for Yios-dropouts, B435, V606, ^775, ^850 are 
considered. 



fled 47 ^850 -dropouts and 27 YiQ 5-dropouts in the HUDF 
main field. IMcLure et al.l (|2012[ ) independently searched 
for galaxies at z > 7 using the photometric redshift tech- 
nique. The objects in their catalog with photometric 
redshifts %>hoto ~ 7 — 9 are well matched with the ob- 
jec ts in the zrso- and yi n5-dropout catalogs constructed 
bv ISchenker etall (|2012h . 

In addition to the dropout galaxies, we stu dy z > 8.5 
star-fo rming galaxy candidates reported by lEllis et aTl 
(j2Q12h . They located all sources by examining the 
stack of the final J125-, «7i4o-? and i^i6o-band images 
and applied the pho tometric redshift technique (see also, 
IMcLure et al.l 120121 ). making use of the full data set ob- 
tained by the UDF12 program and the previous pro- 
grams. They also applied the dropout technique for the 
master catalog, searching for objects undetected at 2a in 
both F105 (> 31.0 mag) and in a combined ACS image. 
By both of the two techniques, they have found seven 
convincing z > 8.5 candidates. 

Morphology measurements for galaxies require a signif- 
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icant detection in not only the central regi on of sources, 
but a lso the outer structures. Recently, M osleh et al.l 
(j2012[ ) reported that, in order to recover the input sizes 
of their realistic simulations, a signal-to- noise (S/N) ra- 
tio of at least 10 is required. To obtain robust esti- 
mates on galaxy morphologies, we set a more strict crite- 
rion for S/N; we analyze our dropouts individually down 
to S/N of 15. The number of zgso-dropouts and Y105- 
dropouts with detection greater than 15a in J125 + J140 
and J140 + Hiqo (about 28.5 mag in the filter- matched 
apertures) is 9 and 6, respectively. 

In order to extend our analysis to fainter magnitudes, 
we divide the samples into three luminosity bins, L 



(0.3 - 1)L* =3 , L 



(0.12 - 0.3)L* =3 , 



L 



(0.048 



0.12)L* =3 , based on their total magnitudes in J125 + J140 
for zgso-dropouts and J140 + #160 fc> r lio5-dropouts. 
Since it is difficult to establish reliable total magnitudes 
for faint sources (S/N < 15) using GALFIT, we subdivide 
the galaxies into luminosity bins based on aperture mag- 
nitudes which contain 70% of a point-source flux, after 
making the appropriate apert ure correction to 10 0% of 
anticipated point-source flux (M cLure et all 120121 ). We 
make median-stacked images separately for the second 
and third brightest luminosity bins. The number of ^850- 
dropouts (Yio5-dropouts) with L = (0.12 — 0.3)L* =3 is 8 
(7), while the number with L = (0.048 - 0.12)L* =3 is 17 
(13). Note that, among the 8 zgso-dropouts (the 7 Y105- 
dropouts) in the second brightest luminosity bin, 7 (3) 
are individually detected at more than 15a in J125 + J140 
( J140 + Hiqo). We do not use stacked images for the ob- 
jects in the brightest luminosity bin, since the numbers 
of the objects are small. We also tried stacking objects 
fainter than L = 0.048L* =3 , but they did not provide 
meaningful size constraints. 
Within the z > 8.5 sample, UDF12-3954-6284 has a 



relatively high photometric redshift, 
the others have gp hoto = 8.6 - 9.5 



nhntn = H.9 while 

(Ell is et all l20ll 



iMcLure et aHl2Q12fi - Thus, we divide them into two sub- 
samples: one with UDF12-3954-6284 and the other with 
the remaining six objects. The average photometric red- 
shift of the latter subsample is about 9.0. Since most of 
these objects are quite faint, we make a stack of the Hiqq 
images of the six z c± 9 objects, giving a detection with 
S/N - 9. 

Note that the nature of the z ~ 12 source is still un- 
certain, beca use of its accompanying diffuse morphology 
(Section 14. 2ft and its luminosity, particularl y given the 
lack of other detections beyond z ~ 10.5 ([Ellis et al.l 
2012). Nevertheless, since no alternative, plausible ex- 
planation for this object has yet been proposed, we ana- 
lyze this object as a z ~ 12 candidate individually. 

In our following analysis, we treat the 9 zgso-dropouts 
and 6 Yios-dropouts with > 15a detections individually, 
and also the 4 stacked objects at z ~ 7 — 8. In addition, 
we analyze the stacked z ~ 9 object and the z = 11.9 
object. 

Note that our bright dropouts have been found in the 
literature as summarized in Table [TJ Table [2] summarizes 
their photometric redshift s reported so far. 



4. SIZES OF GALAXIES AT , 



12 



has the following form: 
E(r) = S e exp | 



l/n 



1 



(1) 



where S e is the surface brightness at the half-light ra- 
dius r e , n is the Sersic index, which is often referred to 
as the concentration parameter; larger n values denote 
steeper inner profiles and highly extended outer wings. 
The parameter r e is the half-light radius, which holds 
half of the total flux inside. To make this definition 
true, the variable b n depends on n. We fit the two- 
dimensional surface brightness profile usin g the GALFIT 
software version 3 (jPeng et alJEOOl l2010l ). which con- 
volves a galaxy model profile image with a PSF profile 
and optimizes the fits using Levenberg-Marquardt algo- 
rithm for x 2 minimization. The output parameters in- 
clude the centroid coordinates of the objects, their total 
magnitude, half-light radius, Sersic index n, axis ratio, 
and position angle. The half-light radius provided by 
GALFIT is the radius along the semi-major axis, a. For 
each galaxy, we calculate the circularized half-light ra- 
dius, r e = ayjbja, where b/a is the axis ratio. The 
initial paramete rs used for profile fittin g are provided 
by SExtractor ([Bertin fc Arnoutsl fl996l ) . and all of the 
parameters, except for the Sersic index, n, are allowed 
to vary during the fitting procedure. The Sersic index 
n is fixed at 1.0, which corresponds to the exponential 
profil^S In this case, b n = 1.678, which is obtained 
by solving the following equation: 7(2n, b n ) = T(2n)/2, 
where 7 is the incomplete gamma function, and V is the 
gamma function. Noise images, required to weight indi- 
vidual pixels in the fit, are taken to be the root mean 
square (rms) maps generated from variance maps pro- 
vided by the data reduction. We also use segmentation 
images which are produced by SExtractor, to mask ob- 
jects other than the object we are interested in during 
the profile fitting. 

4.1. Simulations of Systematic Effects 

Low surface brightness in the outskirts of a galaxy may 
not be correctly measured by GALFIT, leading to system- 
atically low measured half-light radii and/or total mag- 
nitudes. In order to quantify and correct for any such 
systematic effects, we use the following simulations. 

First, we produce galaxy images whose Sersic index 
n is fixed at 1.0, half-light radius r e is randomly cho- 
sen between 0.5 and 10.5 pixels, and total magnitude 
is randomly chosen between 26 and 30 mag. Note that 
axis ratios are fixed at 1 during the simulations. This 
means that the systematic and statistical uncertainties 
will be larger than those obtained by our simulations, 
if output axis ratio is smaller than l P*L Then we con- 
volve them with a PSF image which is a composite 
of bright and unsat urated stellar objects in the HUDF 
([Pirzkal et al.ll2005l ). Figure [1] shows the measured PSFs 
for the J125 + J140, «/uo + #160, and Hi 60 images. The 



The Sersic power law (|Sersidll968f ) is one of the most 
frequently used profiles to study galaxy morphology and 



15 Although Oesch et al. (2 010an set the Sersic index n to be 1.5 
during their analysis, we confirm that measured sizes show very 
little difference if we use n = 1.5. 

16 If axis ratios are parametrized, measurements of circularized 
radii are systematically underestimated by about 10 % and the 
statistical uncertainties in the measured circularized radii are about 
40 % at S/N = 15 (Yuma et al. 2012 in preparation). 
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Figure 1. The contours of PSF images in J125 + «/i40 (left), J140 + H\m (middle), and H±qq (right). The half-light radii of the PSFs are 
0.119 arcsec (3.97 pixels) in J125 + «/i40, 0.124 arcsec (4.14 pixels) in J140 + Hiqq. and 0.123 arcsec (4.12 pixels) in ifi60- 
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Figure 2. Panels (a) — (b), (c) — (d), and (e) — (f) show the results of our simulations for Z850 -dropouts, Yios-dropouts and z > 8.5 
candidates, respectively. These figures show input half-light radius ri in ^ versus output radius ri° ut ^ for a range of output magnitudes, 
26 — 27 (a, c, e) and 27 — 28 mag (b, d, f). The red filled circles and the red error bars denote the average value and the relevant 



(out) 



rms. The blue dashed line shows the relation of r 



(in) 



(out) 



PSF- convolved galaxy images are inserted into empty re- 
gions of the original images before being analyzed in an 
identical manner to the true galaxy sample. 
Figure [2] displays the results of size measurements 

of our simulated galaxies. The panels show vs. 
r (out) each image at two different magnitude ranges 
(26 < ra(° ut ) < 27 and 27 < < 28). We see that 

measurements for all images give low systematic offsets 
for objects with sizes smaller than ~ 4 pixels, although at 
larger sizes the profiles are progressively underestimated 
as the surface brightness of the objects decrease. The 
systematics are also seen to be larger for the fainter ob- 
jects. We also use these simulation results for estimating 



statistical errors in the measurements. 

Figure [3] shows the results for measured total mag- 
nitudes compared to input magnitude. This time the 

results are displayed in two size bins (1 < ri™^ < 3 pix- 
els, 3 < re° ut ^ < 5 pixels) for each image. The results 
for the smaller size bin show that the total measured 
magnitude is robust down to ~ 28 mag. For objects 
fainter than this the measured magnitude is systemati- 
cally fainter than the intrinsic value, and the statistical 
errors increase. The trend is similar for both size bins 
but the results for larger objects show greater systematic 
offsets and statistical uncertainties. 

In summary, our simulations show that GALFIT mea- 
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Figure 3. Panels (a) — (b), (c) — (d), (e) — (f) show the results of our simulations for zsso-dropouts, Yios-dropouts, and z > 8.5 candidates, 

respectively. These figures show input magnitude ra( m ) versus output magnitude m( out ) for a range of output half-light radii, ri ^ =1 — 3 
(a, c, e) and 3 — 5 pixels (b, d, f). The red filled circles and the red error bars denote the average value and the relevant rms. The blue 
dashed line shows the relation of ra( m ) = m( out ) . 



surements of half-light radii and total magnitudes are 
systematically underestimated for faint objects. We cor- 
rect for systematic effects in the half-light radii and total 
magnitudes using the measured offsets in Figures [2] and 
[3j respectively. Note that the errors on r e and total mag- 
nitude reported in this paper are also taken from these 
simulations. 

4.2. GALFIT Measurements 

We perform surface brightness profile fitting for our 
samples at z ~ 7 — 12, using GALFIT and making use 
of our simulation results to correct for any systematic 
effects. We analyze each of the objects with > 15<j de- 
tections individually (9 zgso-dropouts, 6 Yios-dropouts), 
as well as the z = 11.9 object, which is formally detected 
at ~ 8a. We extend the analysis to fainter magnitudes 
using stacked observations. The fainter zg5o- an d Yios- 
dropouts are split into two luminosity bins before stack- 
ing (0.12 < L/L* z = 3 < 0.3 and 0.048 < L/L* z = 3 < 0.12), 
whereas we group all z ~ 9 candidates into a single stack. 

Figure H] presents the results of Sersic profile fitting for 
the 9 bright zgso-dropouts. Shown, from left to right, are 
the 3" x 3" cut-outs of the original image, the best-fit 
model produced by GALFIT, the residual images (origi- 
nal image — best-fit profile) and the segmentation maps 
used for masking all the neighboring objects during the 
profile fitting. Figure [5] similarly shows the results for 
the 6 bright Yios-dropouts. All the objects are cleanly 
subtracted in the residual images. Note, however that 



three of the objects; two of the brightest Z850-dropouts, 
UDF12-4258-6567 and UDF12-3746-6327, as well as one 
of the Yios-dropouts, UDF12-3952-7173 are significantly 
blended with neighboring objects in the original images. 
In addition, one of the Yios-dropouts, UDF12-4470-6442 
shows two cores. The uncertainties in the derived profile 
parameters for these objects will therefore be larger than 
for other isolated objects. 

Figure [6] shows the profile fitting result for the z ~ 12 
object, UDF12-3954-6284. Since the magnitude of the 
z ~ 12 object in i?i6o measured with 0.50"-diameter 
aperture is 29.2 mag, corresponding to S/N ~ 8, the 
profile fitting for this object is quite challenging. Actu- 
ally, the best-fit model galaxy profile seems more elon- 
gated than that in the original image shown in Figure 
[6j which would overestimate of its total magnitude. At 
least, the residual image in Figure [6] does not clearly 
show any noticeable residuals around the central posi- 
tion, although the uncertainties of the fitting parame- 
ters are relatively large as inferred from the moderate 
S/N ratio. If we measure the curve of growth for this 
object, using progressively larger circular apertures, we 
find that the magnitude saturates at 28.8 mag within 
an aperture diameter ~ 0.45". We also find by this ro- 
bust method that the half-light of the source is covered 
by about 0.35^-diameter aperture, and after considering 
the PSF broadening effect, we obtain its half-light ra- 
dius, r*hi = 0.45 kpc, which is consistent with the GALFIT 
measurement within ~ la, and is also nearly equal to 
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Figure 4. Sersic profile fitting results for bright ^850 -dropouts found in the HUDF main field. Shown, from left to right, are the 3 /x x 3 /x 
cut-outs of the original image, the best-fit model profile images, the residual images which are made by subtracting the best-fit images from 
the original ones, and the segmentation maps used for masking all the neighboring objects during the profile fitting. 
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Figure 5. Sersic profile fitting results for bright Y105 -dropouts found in the HUDF main field. Shown, from left to right, are the 3" x 3" 
cut-outs of the original image, the best-fit model profile images, the residual images which are made by subtracting the best-fit images from 
the original ones, and the segmentation maps used for masking all the neighboring objects during the profile fitting. 
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Figure 6. Sersic profile fitting results for the z ~ 12 source, UDF12-3954-6284. Shown, from left to right, are the 3" x 3 /x cut-outs of the 
original image, the best-fit model profile images, the residual images which are made by subtracting the best-fit images from the original 
ones, and the segmentation maps used for masking all the neighboring objects during the profile fitting. 
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Figure 7. 3" x 3" cut-outs of the z ~ 12 source UDF12-3954-6284 from various subsets of the WFC3/IR ffi60-band observations. From 
left to right, the first half of the 2009 dataset, the second half of the 2009 dataset, the full 53-orbit 2009 dataset, the 26-orbit 2012, and 
the full 84-orbit dataset (including 2009, 2012, and other exposures in this field). 
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Figure 8. Same as Figure 2] except that the objects are the stacked Z850-dropouts, whose UV luminosities are L = (0.12 — 0.3)L*_ 3 
(top) and L = (0.048 - 0.12)L* =3 (bottom). 
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Figure 9. Same as Figure |S] except that the objects are the stacked Yios-dropouts. 
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Figure 10. Same as Figure [8] except that the object is the stacked z ~ 9 object. 
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the value reported bv lBouwens et al.l ([2012a[ h ^0.5 kpc. 

Additionally, we note that this object has an unusual 
morphology. It is visually confirmed that the z ~ 12 ob- 
ject has a diffuse filamentary structure stretching from 
north-east to south-west, although the significance is 
very low. This has bee n alread y mentioned very re- 
cently by iBouwens et al.l (|2Q12a| ) . Figure [7] shows the 
cutout H\6o images of this object from various subsets 
and the full data. The diffuse structure is seen in the 
full (2009+2012) data and in the 2009 data. The 2009b 
cutout also shows a low-S/N filament, and the 2009a 
cutout has a similar pattern along the same direction. If 
this diffuse filament is indeed associated with the source 
at z ~ 12, it corresponds to its bright UV continuum 
and/or Lya, which would suggest that this object is ex- 
periencing a major merger event, leading to their high 
star-formation activity. This star formation enhance- 
ment may explain the visibility of such a high-redshift 
galaxy. 

Figure [8] and [9] show profile fitting results for the 
stacked Z850-dropouts and Yios-dropouts, respectively, 
whose UV luminosities are L = (0.12 — 0.3)LJ =3 (top) 
and L = (0.048 - 0.12)L* =3 (bottom). Also shown 
in Figure [10] is the profile fitting result for the stacked 
z c± 9 candidates. Note that we also make averaged (not 
median-stacked) images and perform profile fitting using 
GALFIT, which yield similar fitting results, although for 
some of the average stacks, GALFIT does not provide a 
reasonable fit due to severe confusion with neighboring 
objects. 

In the brightest luminosity bin, L = (0.3 — 1)L*_ 3 , we 
do not perform a stacking analysis, since the numbers 
of the dropouts are small (2 for zgso-dropouts and 3 for 
Yios-dropouts) and the stacked images are significantly 
confused by neighboring objects. Instead, we calculate 
their average sizes and magnitudes; r e = 0.79 ± 0.29 kpc 
and M uv = -20.33 ± 0.20 mag (z ~ 7) and r e = 0.67 ± 
0.28 kpc and M uv = -20.08 ± 0.10 mag (z ~ 8). 

The best-fit parameters are summarized in Table [3] for 
the zgso-dropouts, Table |4] for the Yios-dropouts, and 
Table [5] for the z > 8.5 candidates. The weighted mean 
of half-light radii for the zgso-dropouts and Yios-dropouts 
withL = (0.05-l)L* =3 are 0.35±0.07 kpc and 0.38±0.09 
kpc, respectively. In the next section, we present the size- 
luminosity relation, and investigate the redshift evolution 
of galaxy sizes and SFR surface densities based on these 
results. 

5. RESULTS AND DISCUSSION 

Our measurements of half-light radii in Tables [3]— [5] 
show very small values, typically < 0.5 kpc (see also 
filled symbols in Figure [TT]). The average half-light radii 
of the d ropo uts are only ~ 0.3 — 0.4 kpc at z ~ 7 — 8 
(Section H2j) and at z > 8.5 (Table 0). The half-light 
radius of our z ~ 12 candidate is also remarkably small, 
0.32 ± 0.14 kpc. Even including the la uncertainties, 
these half-light radii are, coincident ally, just as large 
as those of giant molecular associations (GMAs) with 
a mass of ~ 10 7 found in the local universe (e.g., 
Vogel et all 119881 : iRand fc Kulkarnil 119901 : iTosaki et all 
20071 ). 

We investigate the relation between size and luminos- 
ity, i.e. the size-luminosity relation, at each redshift. 
Figure [11] presents the size-luminosity relation for our 
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Figure 11. The size-luminosity relation for 2850 -dropouts (top) 
and Yios-dropouts bottom. The filled circles correspond to the 
stacked objects with UV luminosities of L = (0.12 — 0.3)L* =3 , 
L = (0.048 — 0.12)L*_ 3 , and the averaged values of the objects 
with L = (0.3 — 1)L*_ 3 . The filled squares show the bright objects 
detected at > 15cr without any blending with neighboring sources, 
while the crosses show the objects detected in more than 15a and 
blended with a neighboring source. The dashed curves in each 
figure correspond to a constant star-formation rate density SgFR 
[M© yr- 1 kpc~ 2 ] = 0.1, 0.3, 1, 3, 10 from top to bottom. 

Z85o-dropout and Yios-dropout galaxies at z ~ 7 — 8. 
Our z > 8.5 galaxy candidates are not shown, because 
we cannot constrain the relation with only two measure- 
ments (one from an individual object and one from the 
stack) . In Figure [TTJ fainter galaxies have a smaller half- 
light radius than brighter gala xies. This trend is the 
same as those of local galaxies (|de Jong fc Lacevl |2QQQ[ ) 
as well as high-z dropout gal axies at slightly lowe r red- 
shifts (z ~ 6 — 7), studied by Grazian etaTl (|2012h . Be- 
cause the luminosity of a galaxy depends on two physi- 
cal quantities (surface brightness and size), one needs to 
clarify which quantity is dominant in shaping the size- 
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Table 3 

Surface Brightness Profile Fitting Results for Bright Z850 -dropouts 



Object ID 


RAt 1 
[h:m:s] 


Decl.t 1 
[d:m:s] 


(ap)f 2 

[mag] 


nt 3 


muv 
[ma£ 


f4 


M uv t5 
[mag] 


rJ 6 
[kpc] 






S/N> 15, L/L* =3 


= 0.3 


- 1 


















UDF12-3746-6327 


3:32:37.46 


-27:46:32.7 


27.38 


1.0 


26.25 


± 


0.06 


-20.54 


± 


0.06 


1.09 


± 


0.09 


UDF12-4258-6567 


3:32:42.58 


-27:46:56.7 


27.41 


1.0 


26.74 


± 


0.05 


-20.13 


± 


0.05 


0.50 


± 


0.03 






S/N> 15, 




= 0.12 


-0.3 


















UDF12-4256-7314 


3:32:42.56 


-27:47:31.4 


27.73 


1.0 


27.23 


± 


0.05 


-19.66 


± 


0.05 


0.38 


± 


0.05 


UDF12-4219-6278 


3:32:42.19 


-27:46:27.8 


28.09 


1.0 


27.74 


± 


0.10 


-19.03 


± 


0.10 


0.26 


± 


0.07 


UDF12-3677-7535 


3:32:36.77 


-27:47:53.5 


28.22 


1.0 


27.82 


± 


0.11 


-18.97 


± 


0.11 


0.38 


± 


0.09 


UDF12-4105-7155 


3:32:41.05 


-27:47:15.5 


28.30 


1.0 


27.98 


± 


0.13 


-19.00 


± 


0.13 


0.20 


± 


0.07 


UDF12-3958-6564 


3:32:39.58 


-27:46:56.4 


28.33 


1.0 


27.87 


± 


0.11 


-18.93 


± 


0.11 


0.42 


± 


0.10 


UDF12-3744-6512 


3:32:37.44 


-27:46:51.2 


28.38 


1.0 


28.02 


± 


0.13 


-18.80 


± 


0.13 


0.15 


± 


0.09 


UDF12-3638-7162 


3:32:36.38 


-27:47:16.2 


28.41 


1.0 


28.05 


± 


0.14 


-18.77 


± 


0.14 


0.11 


± 


0.08 



stack 



UDF12z-stackl (L/L* =3 = 0.12 -0.3) 28.27 1.0 27.87 ±0.11 -19.00 ±0.11 0.31 ± 0.09 

UDF12z-stack2 (L/L* z = 3 = 0.048 - 0.12) 29.20 1.0 28.73 ± 0.26 -18.14 ± 0.26 0.36 ± 0.15 

t 1 Coordinates are in J2000. 

t 2 Measured in 0.46^-diameter aperture with the stack of the J125 and J140 images. 
t 3 Sersic index. This is fixed, not measured. 

t 4 Total magnitude measured by GALFIT. The systematic effect is considered. 

t 5 Total absolute magnitude calculated using z p hoto if available, otherwise z = 6.7, which corresponds to the 
peak of the selection function for Z850 -dropouts. 

t 6 Circularized half-light radius r e = a^b/a, where a is radius along the major axis, and b/a is axis ratio. 



Table 4 

Surface Brightness Profile Fitting Results for Bright Y105 -dropouts 



Object ID 


RAt 1 
[h:m:s] 


Decl.t 1 
[d:m:s] 


Op)f2 

[mag] 


nt 3 


rauv 1 " 4 
[mag] 


M uv t5 
[mag] 


r e te 
[kpc] 






S/N> 15, L/L* =3 


= 0.3 


- 1 






UDF12-3879-7071 


3:32:38.80 


-27:47:07.1 


27.21 


1.0 


26.80 ±0.03 


-20.18 ±0.03 


0.36 ±0.02 


UDF12-4470-6442 


3:32:44.70 


-27:46:44.2 


27.69 


1.0 


27.13 ±0.08 


-19.93 ±0.08 


0.60 ± 0.09 


UDF12-3952-7173 


3:32:39.52 


-27:47:17.3 


28.10 


1.0 


26.95 ±0.10 


-20.12 ±0.10 


1.05 ±0.19 






S/N> 15, L/L* =3 = 


= 0.12 


-0.3 






UDF12-4314-6284 


3:32:43.14 


-27:46:28.4 


28.10 


1.0 


27.73 ±0.08 


-19.22 ±0.08 


0.34 ±0.07 


UDF12-3722-8061 


3:32:37.22 


-27:48:06.1 


28.28 


1.0 


27.89 ±0.10 


-19.09 ±0.10 


0.20 ± 0.06 


UDF12-3813-5540 


3:32:38.13 


-27:45:54.0 


28.33 


1.0 


27.99 ±0.11 


-19.23 ±0.11 


0.09 ±0.05 








stack 










UDF12y-stackl (L/L*_ 3 = 0.12 - 
UDF12y-stack2 (L/L* z = 3 = 0.048 


-0.3) 
- 0.12) 


28.40 
29.47 


1.0 
1.0 


28.01 ±0.11 
28.90 ±0.29 


-19.13 ±0.11 
-18.24 ±0.29 


0.35 ±0.13 
0.36 ±0.13 



t 1 Coordinates are in J2000. 

t 2 Measured in 0.50 /; -diameter aperture with the stack of the J140 and if 160 images. 
t 3 Sersic index. This is fixed, not measured. 

t 4 Total magnitude measured by GALFIT. The systematic effect is considered. 

t 5 Total absolute magnitude calculated using z p hoto if available, otherwise z = 8.0, which corresponds to the 
peak of the selection function for Yios-dropouts. 

t 6 Circularized half-light radius r e = a^Jb/a, where a is radius along the major axis, and b/a is axis ratio. 



luminosity relation. We define star-formation rate sur- 
face density, Ssfr, as the average star- formation rate in 
a circle whose radius is r e , 



From equations (|2|) — ((3]), we obtain 



M uv -2.5 log 



Ssfr r\ 



SsFRlMoyr^kpc-^^M. 

irri 



1.4 x 10- 28 • 2 • (lOpc [cm])' 



48.6. 



(2) 



In the case that dust extinction is negligible, a rest-frame 
UV luminosity dens ity approximately correlates with a 
star-formation rate ([Kennicu ttl [l998ah . 



SFR [Mq yr -1 ] = 1.4 x l(r 28 L„ [ergs^Hz -1 ]. (3) 



(4) 

In Figure [TTJ we show constant star- formation rate sur- 
face densities of Ssfr = 0.1, 0.3, 1, 3, and 10 [Mq yr -1 
kpc -2 ] with dashed lines. We find that most of indi- 
vidual galaxies and stacked galaxies fall in the range of 
^sfr — 1 — 10 within their uncertainties. These results 
indicate that both bright and faint z ~ 7 — 8 galax- 
ies have similar star-formation rate surface densities of 
Ssfr — 1 — 10, and that the size-luminosity relation at 
each redshift is mainly determined by the size of galaxies 
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Table 5 

Surface Brightness Profile Fitting Results for z > 8.5 Candidates 



Object ID 


RAt 1 
[h:m:s] 


Decl.t 1 
[d:m:s] 


(ap)f2 
''TV 
[mag] 


nt 3 


muv t4 
[mag] 


Muv 15 
[mag] 


rJ 6 
[kpc] 


UDF12-3954-6284 


3:32:39.54 


-27:46:28.4 


29.24 


1.0 


28.47 ±0.25 


-20.41 ±0.25 


0.32 ±0.14 








stack 










UDF12hz-stack 






29.69 


1.0 


29.11 ±0.49 


-18.21 ±0.49 


0.35 ±0.16 



t 1 Coordinates are in J2000. 

t 2 Measured in 0.50^-diameter aperture with the Hieo image. 

Sersic index. This is fixed, not measured. 
+ 4 Total magnitude measured by GALFIT. The systematic effect is considered. For UDF12-3954-6284, this 
might be overestimated, since the best-fit model galaxy profile seems more elongated than that in the original 
image shown in F igur e [6] If we measure the curve of growth for this source, the magnitude saturates at 28.8 
mag (See Section N~2)) . 

t Total absolute magnitude. We calculate it with z p hoto — 

11.9 for UDF12-3954-6284, considering IGM 
absorption shortward of its Lya wavelength. We use the average photometric redshift, z p hoto — 9.0 f° r the 
stacked object. 

Circularized half-light radius r e = a^b/a, where a is radius along the major axis, and b/a is axis ratio. 



that have a similar star-formation rate surface density. 

We then investigate the size evolution. Since the half- 
light radius depends on luminosity, as displayed by the 
size-luminosity relation, we carefully compare the half- 
light radii of our dropout galaxies within a fixed mag- 
nitude range. Figure [12] presents the average half-light 
radius as a function of redshift for our dropout galaxies 
at z ~ 7-12 with (0.3-l)L* =3 and (0.12-0.3)L* =3 , to- 
gether with dropout galaxies at z ~ 4 — 8 taken from the 
literature. Our measurements of average half-light radii 
are consistent with those from the previous studies at 
z ~ 7 — 8, where we see overlap with previous measure- 
ments. Figure [12] indicates that the average half-light 
radius decreases with redshift from z ~ 4 to 8, which 
is consistent with the claims of previous studies (e.g. 
Ferguson etHI 120041: IBouwens et all 120041 : lHathi et all 



2008; Oesch et al. 2010a). 



UDF12 provides us with the deepest ever near-infrared 
images of the HUDF, allowing our study to extend the 
dynamic range of redshift in the size evolution analysis 
from z~8toz~12, and identifies that the decreasing 
trend holds up to z ~ 12 as shown in Figure [12J if the 
putative z ~ 12 source is real. Although the statistical 
uncertainty of measurement is large, the half-light radius 
of z ~ 12 is r e = 0.32 ± 0.14 kpc in the luminosity bin 
of (0.3 — 1)L* =3 , which is significantly smaller than that 
of z ~ 6 by a factor of 3. Note that we can only plot 
the half-light radius at z ~ 12 on the panel of (0.3 — 
1)L* =3 in Figure [I2j because there are no z ~ 12 galaxies 
with (0.12 - 0.3)L* =3 in the UDF12 data. Similarly, 
our stack of z > 8.5 galaxies have a luminosity fainter 
than (0.12 — 0.3)L* =3 , which is too faint to be compared 
with the baseline of the average half-light radii at z ~ 
4 — 6. However, our results of z > 7 galaxies at these 
faint magnitudes, which are shown as gray filled circles 
in the bottom panel of Figure [12j are consistent with the 
decreasing trend of the half-light radius, albeit with large 
errors. 

This decreasing trend may be explained by the evolu- 
tion of host dark halo radius. According to the analytic 
model in the hierarc hical structure f o rmation framework 
of ACDM (see, e.g., IMo et al . 1998; Mo fc White! l200l 
iFerguson et a l. 2004), the virial radius of a dark matter 



halo is given by 



\100H(zY 



1/3 



(5) 



where H(z) is the Hubble parameter and M v i r is the 
virial mass of the halo. The virial radius is also expressed 
as a function of the circular velocity of dark halo by 



10H(z)' 



(6) 



Since H(z) is approximated by ~ (1 + z) 3 / 2 in a flat 
universe at high redshift s, the redshift evolution of the 
virial radius is r v i r ex H(z)~ 2 / 3 ~ (1 + z)~ x for constant 
halo mass and r v i r oc H(z)~ l ~ (1 + z) -1,5 for constant 
velocity. 

Figure [12] shows the radius-redshift relation of dark 
matter halos for the case of constant halo mass and con- 
stant velocity. Previous studies investigating the radius- 
redshift relation in the red shift range 4 < z < 8 reach 
two different conclusions. IBouwens et al.l (j2004 l2QQ6h 
claim that the relation is roughly (1 + z) _1 , suggestive 
that the sizes of dis k s sca le wit h constant h alo m ass, 
while IFerguson etHI (j2004f ) and lHathi et al.l (|2008f ) ar- 
gue that (1 + z) -1,5 , i.e., the constant velocity case, 
is preferable. We fit the radius-redshift relation over 
a wider range of redshift (extending to z ~ 12) with 
a function of (1 + z) s , where s is a free parameter. 
We take into account our size measurements: for the 
brighter bin, the average sizes of z^q- and Yios-dropouts 
and the size of the z ~ 12 source, and for the fainter 
bin, the measured sizes of the stacks of ^850- and Y105- 
dropouts with L = (0.12 — 0.3)L* =3 . In addition, we 
use the results reported in the litera ture: the average 
sizes at z = 2.5 (Bouwens et al. 2004), the average sizes 
at z = 4 - 6 (jOesch et alJl2010a[ FL We fit the follow- 
ing two functions to the data, logr e = slog(l + z) + a\ 
for L = (0.3 - 1)£*= 3 and logr e = slog(l + z) + a 2 for 
L = (0.12 — 0.3)L* =3 , where 5, a±, and a 2 are free param- 
eters. Varying the three parameters, we search for the 

17 For the fitting , we do not utilize the GALFIT measurements by 
Oesch et al. (2010a), which they did not use as their fiducial ones. 
Note that the fitting result is consistent within lcr, if we include the 
GALFIT measurements instead of their SExtractor measurements. 
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Figure 12. Top: Evolution of the half-light radius across the 
redshift range from z ~ 2 to 12 in (0.3 — 1)L* =3 . The filled circles 
show the average sizes of our zsso-dropouts and Y105 -dropouts, 
and the size of the z ~ 12 object. The open symbols are taken 
from the literature; the open squares a nd triangles are dropout 
gala xies taken from jOesch et al. (2010a), the open diamonds are 
from lBouwens et al.l (|2004T ). After excluding the sample overlaps, 
we fit simple functions of (1 + z) s with the data in both luminosity 
bins, and obtain s = —1.28 ±0.13. which is shown as the solid line. 
The dotted and dashed lines correspond to the case of s = —1.0 
and —1.5, respectively. Bottom: Evolution of the half-light radius 
across the redshift range from z ~ 4 to 8 in (0.12 — 0.3)L* =3 . The 
open and filled black symbols denote the same as those in the top 
panel. The gray filled circles are dropout galaxies in the fainter 
luminosity bin, (0.048 — 0.12)L* =3 . The solid, dotted, and dashed 
lines are the same as those in the top figure. 



best-fitting set of (s, ai, G&2) that minimizes x 2 - The best- 
fit parameters are s = — 1.28±0.13, a\ = 0.99±0.08, and 
d2 = 0.87^0. io- We have checked that exclusion of the 
the putative z ~ 12 object produces no significant change 
in our results, but it is nevertheless interesting that its 
size conforms with the trend established at slightly lower 
redshifts. We note that these results are clearly con- 
sistent with the redshift trend derived by lOesch et al.l 
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Figure 13. Evolution of the SFR surface density Esfr as a 
function of redshift, for dropouts in the brightest luminosity bin, 
L = (0.3 — 1)L*_ 3 . The filled circles correspond to our zsso- 
dropouts and Yios- dropouts. The open squares are taken from 
lOesch et al.l (|201 0a), showing their dropout galaxies, while the 
open circles are the same objects but corrected for dust absorption. 
The dotted line and dashed line shows the case with a constant UV 



luminosity L = L*^ 3 and 0.3L* 



respectively, given that their 



half-light radii follow the simple relati on o f (1 + z) s 
rived in Section \E\ and shown in Figure fT2] 



which is de- 



(|2010at ) from the early UDF09 data, as they reported 
s = —1.12 ± 0.17 for galaxies with luminosities in the 
range L = (0.3 - 1)L* =3 , and s = -1.32 ± 0.52 for the 
fainter galaxies with L = (0.12 — 0.3)L* =3 . However, our 
derived value for s is more accurate, both because of the 
improved data and galaxy samples provided by UDF12, 
and because we have chosen to fit a single value of s to 
both the bright and more modest luminosity galaxies. 

It should be noted again that the above simple constant 
mass or content velocity models assume that the stellar 
to halo size ratio does not change over this redshift range 
(|Mo et al.lll998| ). To properly interpret our result, more 
realistic models are needed which carefully treat the stel- 
lar to halo size ratio, as well as consider effects on galaxy 
sizes from galaxy mergers, torquing, and feedback, based 
on a hierarchical gal axy formation scenario over the full 
redshift range (e.g., ISomerville et al.l [20081 ). These size 
measurements of high-redshift galaxies provide a launch- 
ing point for a theoretical understanding of the structure 
of such galaxies, which has only recently been attempted 
but is of cri tical importance in underst anding their prop- 
erties (e.g.. iMuhoz fc Furlanettoll2012[ ). 

Figure [13] presents the average star-formation surface 
density, Ssfr, as a function of redshift. Note that the 
measurements are shown up to z ~ 8 in Figure [131 be- 
cause the uncertainty of the z ~ 12 measurements are 
too large to place a meaningful constraint. Our galax- 
ies at z > 7 have Ssfr °° 2. Esfr appears to increase 
towards high redshifts. In fact, this increase of Ssfr 
is expected from the decreasing trend of galaxy size at 
a given luminosity (Figure [T2]). Since Ssfr is propor- 
tional to the UV luminosity density in the case of no 
dust extinction, Figure [13] indicates that UV luminosity 
surface brightness is higher for z ~ 7 — 8 galaxies than for 
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z ~ 4 — 5 galaxies by a factor of 2 — 3. Figure [13] has data 
points of dust-corrected Ssfr- The dust-extinction cor- 
rected EgpR is significantly larger than uncorrected Esfr 
at z ~ 4 — 5 ; while a lmost no difference is found at z > 6. 
lOesch et al.l (j2010a[ ) claim that the dust-extinction cor- 
rected Ssfr is roughly constant over the redshift range of 
z ~ 3 — 7. Extending their study, we find this constant 
Ssfr up to z ~ 8. Dotted and dashed lines in Figure 
[13] show the X$fr vames f° r ^=3 and 0.3L* =3 expected 
from the best-fit size evolution of (1 + z) s . These lines 
imply that one would find galaxies with extremely high 
^sfr at z > 10 and beyond. A simple increase of the 
average Xsfr * s n °t expected, however. This is because 
the typical UV luminosity, L*, is fainter at higher red- 
shifts and galaxies with (0.3 — 1.0)L* =3 are quite rare 
at z > 10. In this sense, even higher-redshift galax- 
ies cannot take an extremely high average £sfr beyond 
~ 3 — 10 Mq yr _1 kpc -2 . In the local universe, the 
SFR surface densities of normal disk galaxies are about 
O.OlM© yr _1 kpc -2 , smaller than what we have found for 
z ~ 4 — 8 dropouts. The centers of normal disk galax- 
i es, on the other hand, re ach about IMq yr - 1 kpc -2 
(|Kennicuttlll998bL see also iMomose et al.l[20Toh . which 
is comparable to Esfr of z ~ 4 — 8 dropouts. Because 
local starbursts, especially nuclear st arbursts, have hig h 
Ssfr up to 100-1000M o yr -1 kpc -2 (|Kennicuttlll998b[ ). 
the star-formation surface density of dropout galaxies at 
z ~ 4 — 8 is significantly smaller than those of the ex- 
treme population found in the local universe, which in- 
dicate that star-formation in dropout galaxies is not as 
rapid as that of local extreme starbursts. Speculatively, 
because high-z dropouts are me tal and dust poor galax- 
ies (e.g., iBouwens et al.l l2012b[ ). gas cooling in a given 
amount of molecular clouds of high-z dropouts would be 
less efficient than that of local starbursts. 

6. SUMMARY 

In this paper, we have presented sizes of dropout 
galaxy candidates at z ~ 7 — 12 identified by the 2012 
Hubble Ultra Deep Field campaign. We have stacked 
the new F140W image with the existing F125W image 
and the deeper F160W image, to maximize the available 
depth at rest-frame wavelengths A res t — 1600 — 1700A for 
zgso-dropout and Yios-dropout samples respectively, al- 
lowing secure size measurement from > 15<r detections. 
The extremely deep F105W data ensures that z > 8 
candidates are robust, extending the redshift range of 
reliable objects. We have performed surface brightness 
profile fitting for our samples at z ~ 7 — 12. Our mea- 
surements have shown that the average half-light radii 
of galaxies are very small, 0.3 — 0.4 kpc at z ~ 7 — 12. 
Such sizes are, perhaps coincident ally, comparable to the 
sizes of giant molecular associations in local star-forming 
galaxies. 

Combining our new results at z ~ 7 — 12 with exist- 
ing average size measurements previously reported for 
dropout galaxies at z ~ 4 — 7, we have investigated the 
size evolution of dropout galaxies. We have confirmed 
the trend for size to decrease with increasing redshift 
(at a given luminosity) and have shown that this trend 
appears to extend out to z ~ 12. Motivated by the 
fact that, at least qualitatively, the sizes of the brighter 
(0.3 - 1.0L* =3 ) and somewhat fainter (0.12 - 0.3L* =3 ) 
dropout galaxies show a similar trend with redshift, we 



have attempted to model the size evolution of both sam- 
ples together with a function of the form (1 + z) s over 
the redshift range z ~ 4 — 12. The result is a best- 
fitting parameter value of s = —1.28 ±0.13, approx- 
imately mid-way between the physically expected evo- 
lution for baryons embedded in dark halos of constant 
mass (s = —1) and constant velocity (s = —1.5). This 
evoluti on is consistent with that derived by lOesch et al.l 
(|2010at ). albeit our derived evolution with redshift is 
slightly steeper than that derived bv lOesch et al.l (|2010af ) 
for the brighter galaxies. We have checked that our best- 
fitting value of s is not significantly affected if the puta- 
tive z c± 12 galaxy is excluded, but it is interesting that 
this object has a half-light radius which is perfectly con- 
sistent with our best-fitting relation. 

We have also found that a clear size-luminosity rela- 
tion, such as that found at lower redshift, is also evident 
in both our zs5o~ and Yios-dropout samples. This rela- 
tion can be interpreted in terms of a constant surface 
density of star formation over a range in luminosity of 
0.05 — 1.0L* =3 . These size-redshift and size-luminosity 
relations suggest that galaxy sizes at z > 4 are not simply 
decided by the evolution of the constant mass or veloc- 
ity of the parent halo and/or follow in the evolution of 
the stellar to halo size ratio with a similar star-formation 
rate density. 

Finally, our results also strengthen previous claims 
that the star-formation surface density in dropout galax- 
ies is broadly unchanged from z c± 4 to z ~ 8 at 
XIsfr — 2 M yr -1 kpc -2 . This value is 2 — 3 orders 
of magnitude lower than that found in extreme starburst 
galaxies, but is very comparable to that seen today in 
the centers of normal disk galaxies. This provides fur- 
ther support for a steady smooth build-up of the stellar 
populations in galaxies in the young universe. 
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